home *** CD-ROM | disk | FTP | other *** search
/ Languguage OS 2 / Languguage OS II Version 10-94 (Knowledge Media)(1994).ISO / gnu / glibc108.zip / glibc108 / sysdeps / m88k / m88110 / mul_1.S
Text File  |  1994-01-15  |  2KB  |  81 lines
  1. ; mc88110 __mpn_mul_1 -- Multiply a limb vector with a single limb and
  2. ; store the product in a second limb vector.
  3.  
  4. ; Copyright (C) 1992, 1994 Free Software Foundation, Inc.
  5.  
  6. ; This file is part of the GNU MP Library.
  7.  
  8. ; The GNU MP Library is free software; you can redistribute it and/or modify
  9. ; it under the terms of the GNU General Public License as published by
  10. ; the Free Software Foundation; either version 2, or (at your option)
  11. ; any later version.
  12.  
  13. ; The GNU MP Library is distributed in the hope that it will be useful,
  14. ; but WITHOUT ANY WARRANTY; without even the implied warranty of
  15. ; MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  16. ; GNU General Public License for more details.
  17.  
  18. ; You should have received a copy of the GNU General Public License
  19. ; along with the GNU MP Library; see the file COPYING.  If not, write to
  20. ; the Free Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA.
  21.  
  22.  
  23. ; INPUT PARAMETERS
  24. ; res_ptr    r2
  25. ; s1_ptr    r3
  26. ; size        r4
  27. ; s2_limb    r5
  28.  
  29. #include "sysdep.h"
  30.  
  31. ENTRY (__mpn_mul_1)
  32.     ld     r6,r3,0
  33.     sub     r4,r0,r4
  34.     sub     r3,r3,r4        ; r3 is offset s1_ptr
  35.     sub     r2,r2,r4
  36.     sub     r8,r2,8        ; r8 is offset res_ptr
  37.     mulu.d     r10,r6,r5
  38.  
  39.     addu     r4,r4,1
  40.     bcnd     eq0,r4,Lend
  41.      addu.co r2,r0,0        ; clear cy_limb
  42.  
  43. Loop:    ld     r6,r3[r4]
  44.     addu.cio r9,r11,r2
  45.     or     r2,r10,r0        ; could be avoided if unrolled
  46.     addu     r4,r4,1
  47.     mulu.d     r10,r6,r5
  48.     bcnd     ne0,r4,Loop
  49.      st     r9,r8[r4]
  50.  
  51. Lend:    addu.cio r9,r11,r2
  52.     st     r9,r8,4
  53.     jmp.n     r1
  54.      addu.ci r2,r10,r0
  55.  
  56. ; This is the Right Way to do this on '110.  4 cycles / 64-bit limb.
  57. ;    ld.d    r10,
  58. ;    mulu.d
  59. ;    addu.cio
  60. ;    addu.cio
  61. ;    st.d
  62. ;    mulu.d    ,r11,r5
  63. ;    ld.d    r12,
  64. ;    mulu.d    ,r10,r5
  65. ;    addu.cio
  66. ;    addu.cio
  67. ;    st.d
  68. ;    mulu.d
  69. ;    ld.d    r10,
  70. ;    mulu.d
  71. ;    addu.cio
  72. ;    addu.cio
  73. ;    st.d
  74. ;    mulu.d
  75. ;    ld.d    r10,
  76. ;    mulu.d
  77. ;    addu.cio
  78. ;    addu.cio
  79. ;    st.d
  80. ;    mulu.d
  81.